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Abstract 

In this paper, we perform a detailed study of the scaling properties of a ferromagnetic thin film model. 
Recently, interest has increased in the scaling properties of the magnetic domain wall (MDW) motion in 
disordered media when an external driving field is present. We consider a (l+l)-dimensional model, based 
on evolution rules, able to describe the MDW avalanches. The global interface width of this model shows 
Family-Vicsek scaling with roughness exponent ( ~ 1.585 and growth exponent /3 cs 0.975. In contrast, this 
model shows scaling anomalies in the interface local properties characteristic of other systems with depinning 
transition of the MDW, e.g. quenched Edwards-Wilkinson (QEW) equation and random-field Ising model 
(REIM) with driving. We show that, at the depinning transition, the saturated average velocity Vsat 
vanished very slowly (with 0 ~ 0.037) when the reduced force / = p/pc —1 —^ O’*'. The simulation results show 
that this model verifies all accepted scaling relations which relate the global exponents and the correlation 
length (or time) exponents, valid in systems with depinning transition. Using the interface tilting method, 
we show that the model, close to the depinning transition, exhibits a nonlinearity similar to the one included 
in the Kardar-Parisi-Zhang (KPZ) equation. The nonlinear coefficient A ^ with (j) ~ —1.118, which 
implies that A ^ 0 as the depinning transition is approached, a similar qualitatively behaviour to the driven 
REIM. We conclude this work by discussing the main features of the model and the prospects opened by it. 

1 Introduction 

The jerky motion of magnetic domain walls (MDW) in ferromagnetic media, when an external driving field 
is present, is basically the resultant of magnetic interactions, the internal disorder of the medium and the 
thermal fluctuations. The quenched disorder (or Barkhausen noise) in the media is a consequence of vacancies, 
defects, impurities and dislocations, among others. The interaction between the molecular magnetic field 
and the magnetic dipoles locally governs the movement of the MDW between multi-metastable states of 
quiescence. The dominant dipole-dipole interactions -dipolar (or long-range) and exchange (or short-range)- 
and the thermal fluctuations favour the elastic displacements of the wall, which contributes to moving it from 
the quiescence state. The jerky movement of the MDW between quiescence states is known as Barkhausen 
avalanches or jumps [I]. This complex microscopic dynamic leads to unexpected properties observed both 
experimentally and theoretically. 

The study of ferromagnetic films with Barkhausen noise has reopened the study of general systems with 
avalanches and the discussion about universality classes currently accepted. Such systems can be classified 
according to the values of the scaling exponents obtained from experiments and models. Systems with mag¬ 
netic avalanches are classified by measuring the power-law exponents [5] of the Barkhausen distributions of 
avalanche-size r and -duration a, and the average of avalanche-size as a function of the avalanche-duration 
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Ifavz, and/or the power spectrum d. The studies have led us to accept that there are two universal¬ 
ity classes PHS] associated to the range of the dominant interactions (dipolar or exchange), which prevail in 
amorphous and polycrystalline materials, respectively. For materials that exhibit three-dimensional magnetic 
behaviour (bulk materials including ribbons and sheets) these two universality classes have been accepted, 
even when the experimental values of the exponents were widely dispersed within each class. From the begin¬ 
ning, the experimental study of ferromagnetic films has shown dissimilar results according to the techniques 
employed. The magneto-optical techniques have led to the conclusion that when decreasing the film thickness, 
the two universality classes are maintained even though the range of exponent values change significantly. 
The results based on magneto-optical measurements (restricted to distributions of avalanche-sizes), cross- 
referenced to theoretical predictions, confirm a two-dimensional magnetic behaviour of the thin films [ZHIe]; 
these results are independent of the range of film thickness. This conclusion is consistent with the fully ac¬ 
cepted assertion that systems which exhibit universality are independent of the dynamic microscopic details 
and are controlled by dimensionality and range of interactions, among others properties. Recent reports, 
based on measurements via the inductive technique, suggest that the two-dimensional behaviour accepted 
for films cannot be generalized for all thickness ranges [nno]. Some of these experimental results for the 
scaling exponents grouped films with different thickness in the three-dimensional universality class. In the 
light of this existing disagreement, we think it is necessary to add to the traditional view another approach 
followed in the study of interfaces in disordered media. Thus, for ferromagnetic systems with avalanche 
dynamics there will be at least two criteria for the classification of universality. We focus on the study of 
the scaling properties of the MDW in a disordered medium through a ferromagnetic thin film model. We 
establish the exponents characterizing the pinning-depinning transition and their scaling relations. 

Driven interfaces in disordered media at the depinning transition can be described by stochastic equations 
of motion or lattice growth models. Among the local equations, the quenched Kardar-Parisi-Zhang equation 
(QKPZ) describes systems that can be grouped into two universality classes for their properties at the 
depinning transition; isotropic and anisotropic. The two classes differ fundamentally in the coefficient 
behaviour of the nonlinearity near the transition, as shown below. Some lattice models that can be identified 
with one of these two classes of universality, e.g. random field Ising model with driving belongs to the 
isotropic universality class while directed percolation depining (DPD) models belong to the anisotropic 

universality class [531[21]. There are other models which cannot be included in a QKPZ description type, 
such as the interface evolution models in fractal media, which are grouped in the called isotropic percolation 
depinning (IPD) universality class [3^. The QKPZ equation for the height h = h{x,t) is 

^ = F + i^V^h++ T]{x,h) , (1) 

where F is an external driving force and rj is the quenched noise with mean value equal to zero and correlation 
{r]{x,h)ri{x',h')) = 2D5{x — x') A{h — h'). Here A is an even function which decreases monotonically, 
decaying rapidly to zero at a characteristic distance a, and A = d if a = 0. In the anisotropic universality 
class at the depinning transition the nonlinear coefficient A diverges when the reduced force / —>■ O’*". In 
contrast, the isotropic class A —?► 0 when the / ^ O'*'. Taking A = 0 in equation o, the quenched Edwards- 
Wilkinson (QEW) equation is achieved, and the isotropic class is usually called QEW universality class. 

There are two possible theoretical approaches to explain the depinning transition of the driven MDW 
in a disordered medium: the QEW equation |26H29) and the microscopic lattice models with interactions 
and disorder, e.g. driven RFIM [5DH55] . The QEW equation is a phenomenological coarse-grained equation 
which does not take into account interactions and structure. The MDW described by the QEW equation 
is a single-valued elastic surface plUtdti) . In contrast, the MDW described by the known lattice models is 
not single-valued as a result of islands and/or overhangs [37l|38]. Equations and models can be classified 
according to their scaling properties (usual and anomalous), taking into account the relationship between 
the global, local and spectral roughness exponents [5^MT] . The QEW equation and driven RFIM belong to 
the anomalous scaling class. The subclass of anomalous scaling (intrinsic, super-rough, faceted, or other) to 
which these systems belong is currently still under study [20HS21I1M1H1. 

In a recent paper Muraca and Buceta (BM) have introduced an evolution model of the MDW that takes 
into account only short-range interactions, in order to describe avalanches in a ferromagnetic thin film with 
disorder when a driving field is present [0]. They have found that at the depinning transition the model 
predicts r = 1.29 ± 0.02 and a = 1.55 ± 0.05, which are similar to the results obtained experimentally in 
hlms where exchange interactions dominate [numnnD]. In our previous work m we showed that the 
BM model has anomalous scaling in its local statistical observables, similarly to the QEW equation and the 
driven RFIM but with different anomalous class. The BM model has intrinsic anomalous scaling, while the 
QEW equation has super-rough anomalous scaling and the driven REIM does not have an anomalous 
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scaling univocally defined [12]. Given these interesting properties we decided to study the global scaling at 
the depinning transition using a simplified version of the BM model introduced by us, which maintains all 
the properties listed above. In Section 2, we present the lattice model under study. In Section 3, we present 
the theoretical framework that supports the results obtained from the Monte Carlo simulations performed, 
and establish scaling relations between the exponents of the global roughness and average velocity associated 
with the depinning transition. In Section 3, we also analyse the simulation results and show how these results 
are in agreement with the theory. In Section 4, we show how our model, at the depinning transition, behaves 
similarly to other systems that are within the isotropic universality class. Using the method of the tilted 
surface, the simulation results show that near the transition there is a QKPZ nonlinearity, which vanishes 
at criticality. Finally, a summary is given in Section 5. 


2 The model 

We consider a portion of the ferromagnetic thin-film that includes only two magnetic domains separated by 
a MDW. On either side of the wall we assume opposite macroscopic magnetizations in the easy direction. 
We suppose that the medium is composed of magnetic dipoles and point defects randomly distributed, both 
arranged in the nodes of a 2-dimensional square lattice of edge length L. The BM model assumes that 
two defects cannot be first-neighbours to each other. In contrast, in this model we assume that defects 
are arranged in a completely random way without restrictions. The wall in its movement does not include 
defects. The structure of the MDW is considered merely as a monolayer of dipoles perpendicular to the 
easy direction. The model takes into account only exchange interactions between nearest-neighbour (NN) 
dipoles. To simulate the movement of the MDW, we assume that the lattice has periodic boundary conditions 
and that point defects can be represented by a random pinning force ??(*, j), uniformly distributed in [0,1], 
assigned to each node (i,j) of the lattice, with i = 1,... ,L and j G N. Taking a lattice with density p of 
dipoles and I — p of defects, if < p the node {i,j) has a dipole; otherwise, it has a point defect. We 

characterize the disorder in the lattice using the function F{i,j) = 0(p — r]{i,j)), where Q{x) = I if a: > 0 
and 0(x) = 0 if X < 0. If {k,nk) is a node of the MDW with k = 1,..., T, as the wall is a monolayer of 
dipoles F{k,nk) = 1- Otherwise, any dipole (or defect) outside the MDW is described by F{k,£) = 1 (orO) 
with i ^ rik- 

The evolution rules of this model consider only ferromagnetic exchange (or short range) interactions, 
taking into account the balance of the magnetic moment on either side of the MDW in a neighbourhood of 
the point which is evolving. When there is a local unbalance in the opposite direction of the movement, the 
wall searches the equilibrium with probability c. However, if there is balance, with probability 1 — c, the 
wall can only move if there is external forcing which can remove it from its metastable state. We start the 
Monte Carlo simulation with a flat wall, i.e. initial condition rij = 1 for all i. We introduce, at a fixed time 
on the MDW, the relative position of neighbouring nodes respect to the node of the selected column: 

Xj = rij+i - Uj , yj = - Uj . ( 2 ) 

At time t any column j is chosen randomly. The MDW will advance from the selected node (j, rij ) to the 
node (j, n'j), in the interval St, according to simple evolution rules. Those rules are frustrated if F{j, n') = 0 
or F{j, n'j — 1) = 0; in the first case because the MDW is composed only by dipoles, and in the second case 
because the MDW needs a larger unbalance of energy to overcome a defect. The evolution rules are 

1. With probability c, if Xj +yj > 2 the selected node is moved one unit (i.e. n'j = Uj + 1); otherwise the 
selected node is pinned (i.e. n'j = Uj). 

2. With probability I — c, if Xj = yj > 0 the selected node is moved one unit above its NN columns 
(i.e. n'j = njj .1 -|- I); otherwise the selected node moves (or not) to the maximum between the selected 
column and its NN columns (i.e. n'j = max(nj_i,rij+i)). 

These rules applied over a lattice with randomly distributed defects lead to the same properties as those 
obtained by the BM model: avalanches and scaling anomalies, with similar quantitative conclusions [6ll4^. 

3 Scaling properties 

The scaling theory applied to driven interfaces in disordered media has developed strongly in the last three 
decades, starting from models and equations to describe the phenomenon of pinning-depinning transition |51j . 
Such phenomenon is observed in the MDW motion |52] , the imbibition of fluids in porous media [53j and flux 
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p 

9 

t'h 

/3 

/3c 

0.88700 

0.0415(5) 

1.6206(151) 

1.0228(11) 

0.9750(5) 

0.88710 

0.0389(6) 

1.5224(116) 

0.1005(12) 

0.9751(5) 

0.88715 

0.0376(6) 

1.4730(100) 

0.9777(1) 

0.9751(5) 

0.88716 

0.0374(6) 

1.4631(97) 

0.9749(2) 

0.9751(5) 

0.88717 

0.0371(6) 

1.4532(94) 

0.9699(4) 

0.9751(5) 

0.88720 

0.0364(6) 

1.4234(85) 

0.9565(5) 

0.9751(6) 

0.88730 

0.0338(7) 

1.3233(65) 

0.9195(12) 

0.9751(6) 


Table 1: We show the measured values of the power-law exponents (columns 2-4) as a function of the 
control parameter p* very close to the depinning transition for a lattice size L = 10^. Columns 2 and 3 
show the exponents 9 and t'h of the Usat and Wsat, respectively, as a function of reduced force / = p/pc — 1, 
with Pc — p*, measured in the interval [10"^, 3 x 10^]. Column 4 shows the exponent /3 of the IT as a 
function of time t measured in the same interval. Column 5 shows the growth exponent /3c calculated by 
equations ([5]) and (jH]) with the exponent data 9 and (columns 2-3). Notice that the best approximation 
to the critical value Pc is obtained for L = 10“^ with p* = 0.88716. 
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Figure 1: (Color online) Left: log-log plot of the average saturation velocity Usat as a function of the 
reduced force / = p/pc — 1 close to the depinning transition, with pc — p* = 0.88716. The measured 
exponent is 9 = 0.0376(6) (see Table [T]). Right: log-log plot of the saturation global interface width 
Wsat as a function of the reduced force / close to the depinning transition. The measured exponent is 
= 1.4631(97) (see Tabled]). We take the parameters p > Pc and c = 0.1, and lattice size L = 10"^ in 
the simulations corresponding to both graphs. 
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motion in superconductors [M], among others. The theoretical studies are based on the average velocity of 
the interface, which plays the role of order parameter of the pinning-depinning transition. We assume that 
the system enables a control parameter p, which takes a critical value Pc^ such that if p < Pc the interface 
becomes pinned by disorder and, otherwise it is able to move indefinitely with an average velocity u(/), 
where / = p/pc — 1 is the reduced force. When p p+, the average saturation velocity Usat vanishes as 

(3) 

where 0 > 0. Near the criticality Pc, large and finite regions of the interface are pinned by disorder. At the 
transition, there are correlation lengths (and correlation time) of these pinned regions that diverge as 

\fr^, ( 4 ) 

where t'j > 0 is a correlation length (or time) exponent. Here refers to the lateral (or parallel) correlation 
length ^x, the vertical (or perpendicular) correlation length or the correlation time At the criticality 
p = Pc, in a finite system with lateral size L, the dynamic scaling hypothesis is that the global interface 
width W verifies the Family-Vicsek scaling relation m 

W{L,t)^L^ F{t/L^), (5) 

where C and z are the roughness and dynamic exponents, respectively. The scaling function F{x) ^ when 
a: ^ 1 and goes to a constant when a: ;§> 1, where the growth exponent /? fulfils the relation 

C = Pz. ( 6 ) 

Also, at the criticality for a large and finite system, the average velocity verifies the scaling relation 

u(L,t)~T-^^G(t/L^), (7) 

where 6 is the velocity exponent. The scaling function G{x) ~ x~^ when a; <C 1 and goes to zero when 
a; ^ 1. The exponents that characterize the pinning-depinning transition [by eqs. ([2|) and (|H)] are related to 
the exponents 6 , C and z in a simple way, as we show below. 

Near the depinning transition, large regions of the interface are pinning. The full pinning is achieved 
when the lateral correlation length ^x is of the order of the system size L (i.e. ^x ~ A), the vertical correlation 
length is of the order of the global interface width W (i.e. bF) and the correlation time is of the 
order of the crossover time tx (i.e. G)- Here, the crossover time G is the time at which the correlations 

are propagated throughout the system, saturating the global interface width, thus tx ^ L^. Eq. 0 shows 
that ^ which allows us to obtain that W ^ Taking into account that W L'^ when t tx 

[see eq. 0], we conclude that t'h = C ^x- Also eq. © shows that ~ which allows us to obtain 

that W ~ Considering that W t^ when t ~ tx we conclude that = f] vt. Thus, we obtain the 

following relations between correlation exponents 


Vh = Cvx = 13 vt . ( 8 ) 

Close to the criticality Pc but below it, the average velocity of the interface v ^ t~^ when t ~ tx, and later 
the velocity goes to zero quickly [see eq. ©]• In terms of the reduced force / the velocity v |/|'^'^‘. Taking 
into account eq. © we obtain the relation 9 = Svt if / > 0 m- Assuming that the stationary average 
velocity Vst ■Ch/'Ct E is easy to conclude that Vst ~ and 6 = vt — if / > 0. Thus, we obtain the 

following additional relations between correlation exponents 

9 = i^t - Vh = S i^t . (9) 

Intensive Monte Carlo simulations of our model with large systems close to the depinning transition indicate 
that the growth exponent is /3 ^ 1, which is confirmed through the calculus of the power-law exponent of 
the global interface width W and the saturation height-height correlation function G2,sat [IS]- To set /3 
as accurately as possible it is first necessary to establish the value of p which is closest to its value at the 
depinning transition. A usual procedure is to determine the critical value Pc by epidemic analysis of 

the velocity as a function of time. At criticality ii t the velocity v ^ t~^. From equations © and © 

we can deduce that velocity and growth exponents are related by 

S + P = l, (10) 
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Figure 2: (Color online) Left: plot of as a function of time t. We show the simulation results for 

c = 0.10 and L = 10"*’ and several values of p very close to critical value. Notice that for p = 0.88716, very 
close to the depinning transition, the curve approaches a constant in the interval [10^,3 x 10^]. Right: 
log-log plot of IF as a function of time t. The double arrow o shows the interval where the growth 
exponent (3 = 0.9749(2) is measured. 



Figure 3: (Color online) Log-log plot of the velocity v as a function of time t with lattice size L = 10"* and 
c = 0.10, for several values of p arround the critical value. When p = 0.88716 (red diamond symbols) very 
close to criticality we observe the power law behaviour v ^ with S = 0.0113(62) which is measured 
in the interval [10^, 3 x 10^]. 


equation that is strictly valid in the thermodynamic limit at the depinning transition. By accepting this 
equation we deduce that ^ ^ 0 for our system, which makes the application of the epidemic analysis method 
a disadvantage. Alternatively, the suitable power-law behaviour of the saturation velocity Usat and the 
saturated interface width Wsat as a function of the reduced force / provides a method to establish the 
critical value Pc- Table [1] shows the measured exponents 6 of Usat and Vh of ITsat for several values of the 
control parameter p* close to the depinning transition. With these data we calculate the growth exponent 
/3c — 0.9751 for different values of p*. In contrast, the direct measurement of /3 shows that the value /3c 
is reached for p* = 0.88716 very close to the depinning transition, with /3 ~ 0.9749. Figure [1] shows the 
power-law behaviour of Usat and Wsat as a function of / very close to the critical value Pc- Figure [2] shows 
the interval where the critical power-law behaviour of IF as a function of time t is measured. Figure [3] shows 
the power-law behaviour of the velocity as a function of time, with measured exponent 6 = 0.0113(62) for 
L = 10^, close to the transition. With this method, we find S + ^3 = 0.9862(64) for the same lattice size, a 
very good approximation to the equation (flUl) . We determine C = 1.5851(68) from the saturation interface 
width IFsat as a function of size L for p = 0.88716 (see Figure 0]). With the measured values from Tabled] 
we obtain, by equation ([8]), vt = 1.5008(103) and = 0.9230(100). 

The fact that /3 ^ 1 (or 0^0) has important consequences. The validity of equation (|3|) is confirmed 
because 9 is never zero. Thus, the interface velocity vanishes very slowly when the system approaches 
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Figure 4: (Color online) Log-log plot of the saturation global interface width 114at as a function of lattice 
size L with p = 0.88716 close to criticality and c = 0.10. We measure the power law behaviour Wgat ~ 
with C = 1.5851(68). 


criticality. This hallmark places our model in the limit of the known systems with depinning transition. 
Notice that as \ then, at the transition, the vertical correlation length and the correlation time 

diverge similarly. In contrast, the DPD models verify that vt ~ i/y. [from eq. ([5]) with C ~ /3]; then, at the 
transition, the lateral correlation length and the correlation time diverge similarly. 


4 Behaviour at the depinning transition 

The method of measuring the saturation velocity of the tilted interface, for different values of the reduced 
force /, is useful not only to determine QKPZ nonlinearities but also to establish the universality class at the 
depinning transition. For a fixed value of / when the interface velocity depends on the slope s the system 
includes nonlinearities. The coarse-grain limit only maintains quadratic nonlinearities, which are evident 
with slopes s <S 1. The average saturation velocity of the tilted interface is 

VsAs, f) = Vse^tif) + Hf) s'^ , ( 11 ) 

where Vsat{f) is the saturation velocity of the untilted surface [eq. (|31)] and A(/) is the coarse-grained nonlinear 
coefficient. The velocity of driven interfaces in disordered media, near the depinning transition (/ O"*"), 

shows two opposite behaviours. In systems included within the anisotropic universality class A diverges when 
/ goes to zero. In contrast, other systems are included within the isotropic universality class at the depinnig 
transition when A approaches zero. Simulations of several models show that 

A^/-^. ( 12 ) 

At the depinning transition, if the tilt exponent (j) > 0 {(f) < 0) then the coefficient A diverges (vanishes) and 
the model belongs to the anisotropic (isotropic) universality class. Substituting equations Q and (fT^ into 
equation (fTTl) we find [59] 

Vsat{s,f)(xf+-/f~^s'^, (13) 

with 7 > 0 a constant and 9 > 0. Numerical simulation results of different systems show that Usat is an 
increasing function of /, i.e. dusat/d/ oc [0 — ((>7 > 0. In systems included within the 

anisotropic universality class ^ > 0 and 0 ~ [MSS; then, Vsat always increases if = (^/7 

where Sc is the cutoff slope. Instead, in systems included within the isotropic universality class (at the 
depinning transition) ^ < 0, then Vsat always increases with /. Introducing the crossover slope Sx such 
that Sx ~ the systems included in these universality classes show the following scaling for the tilted 

saturation velocity 

/ „2 \ 

Vsat{sJ)-^fvi^j^j, (14) 
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Figure 5: (Color online) Left: Plot of Usat as a function of the tilt s for several values of p above the 
depinning transition. We note quadratic dependence for s < 1. Right: log-log plot of the scaling function 
y as a function of x [see equation (fHl) ] for several values of / = p/pc — 1, with pc — p* = 0.88716. We 
take 9 = 0.0376 and 6 + (j) = —1.0801. 


where scaling function V{x) ~ x when x 3> 1 and goes to a constant when x <C 1. Introducing the 
characteristic slope s = ^h/Cx it is possible to relate the exponent (j) to other known exponents [23]. If s ^ Sx 
[as ^ and ^ then, for / > 0 we conclude that 

9 + (j) = 2{v^ - z/h) . (15) 

Substituting equation (I5|) into equation (fT51) we obtain 

0 + <(> = 2i/hQ-l) . (16) 

With the measured values of z^h, 9 and C for pc — p* = 0.88716 (from Table[T|) we obtain cj) ~ —1.1177(157) 
[or 9 + (j) = —1.0801(151)]. We use these values (see Figure[5|) to show that the scaling relation [equation (fHl) ] 
for the tilted saturation velocity is verified in our model. The negative value of (j) means that A —> 0 when 
/ —>■ 0+ [from equation (fl^ j. a characteristic behaviour of driven RFIM and QEW equation. 


5 Conclusions 

In summary, in this work we perform a detailed study of the scaling properties of the driven MDW in a 
disordered medium at the depinning transition for a ferromagnetic thin film model. We have chosen this 
model for the study of scaling since it has proven to be highly predictive. First, the model has demonstrated 
to be in agreement with the experimental results of the distributions of avalanche-size and -duration. Second, 
this model has predicted scaling anomalies such as the driven RFIM model and QEW equation. In this paper, 
we apply the known scaling theory that has been developed to understand the motion of driven interfaces in 
disordered media. We have found that the driven MDW verifies the scaling relations established, although 
with some peculiarities. We show that at the depinning transition (p = pc) the relation between scaling 
exponents P + S = 1 holds. Very close to the depinning transition, we show that the growth exponent 
/3 ^ 1, which is a particular feature of this model. As a consequence, the vertical correlation length and 
the correlation time diverge similarly at the depinning transition. Furthermore, the average saturation 
velocity as a function of time shows a power-law behaviour that decays very slowly, approaching a constant, 
since the velocity exponent 5 ~ 0. Also, the average saturation velocity of the interface vanishes very slowly 
when the reduced force goes to zero due to the exponent 0^0. The values of these exponents differ from 
the range of values obtained by other authors for driven RFIM and QEW equation [^[^. Still, we show 
that this model verifies all relations between scaling exponents that are valid for systems with depinning 
transition. We conclude that our model, at the depinning transition, behaves similarly to systems that are 
within the isotropic universality class. Results of simulating our model with the tilted surface show that there 
is a KPZ nonlinearity, which vanishes at criticality. The tilt exponent determined by the simulations confirms 
the known relation with other scaling exponents. When our model approaches the depinning transition we 
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observe a similar qualitatively behaviour to the driven RFIM. However, the values of the exponents of the 
two models differ appreciably. In contrast, the values of the exponents approach each other in the driven 
RFIM and the QEW equation. We think that this difference in the values of the exponents of our model 
with others could be given by correlations of quenched disorder in the motion of a MDW. This conjecture is 
suitable taking into account the results of driven RFIM with different forms and strengths of disorder, which 
lead to a wide dispersion of the exponent values [30] ■ Based on these conclusions, future studies should focus 
on a more general classification of universality to describe the statistical properties of the MDW moving in 
disordered media. We believe that this discrete model is a good candidate to be studied further in order to 
obtain the Langevin equation that characterizes the universality class of the system, by using coarse-grained 
techniques [64ll65] . 
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